Anomalous Nernst effect from a chiral d-density wave state in underdoped cuprate superconductors 
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We show that the breakdown of time-reversal invariance, confirmed by the recent polar Kerr effect measure- 
ments in the cuprates, implies the existence of an anomalous Nernst effect in the pseudogap phase of underdoped 
cuprate superconductors. Modeling the time-reversal-breaking ordered state by the chiral d-density-wave state, 
we find that the magnitude of the Nernst effect can be sizable even at temperatures much higher than the su- 
perconducting transition temperature. These results imply that the experimentally found Nernst effect at the 
pseudogap temperatures may be due to the chiral d-density wave ordered state with broken time-reversal invari- 
ance. 
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I. INTRODUCTION 

Even after two decades of intensive research, the physics 
of the high temperature cuprate superconductors is as elusive 
as ever—. The principal mystery surrounds the underdoped 
regime, which evinces a well-formed quasiparticle gap even at 
temperatures well above the superconducting transition tem- 
perature, T c . The recent observation of a non-zero polar 
Kerr effect (PKE) in the underdoped YBCO 2 , which demon- 
strates macroscopic time-reversal (TR) symmetry breaking in 
the pseudogap phase, is a step forward in solving the pseu- 
dogap puzzle. The PKE appears roughly at the same tem- 
perature, T*, where the pseudogap develops 2 . Near optimum 
doping, the PKE appears at a temperature below T c , consis- 
tent with the existence of a zero temperature quantum phase 
transition under the superconducting dome. This observation 
suggests that the TR symmetry breaking and the pseudogap 
in the cuprates may have the same physical origin, which is 
also unrelated to the d-wave superconductivity itself. Sim- 
ilar conclusion was also reached eariler by muon spin rota- 
tion experiments 3 . In this work we predict the existence of 
an anomalous Nernst effect associated with the TR symme- 
try breaking which should be present along with the observed 
PKE in the underdoped cuprates. Our results demonstrate that 
the existence of anomalous Nernst effect at temperatures as 
high as the pseudogap temperatures (see below), where the 
vortex excitations of the superconductor are unlikely to be 
present, may imply an ordered state with broken TR symme- 
try in the pseudogap regime of the underdoped cuprates. 

It was proposed earlier 4,5 that the id x 2_ y 2 density-wave 
(DDW) state may be responsible for the pseudogap behavior 
in the underdoped cuprates. In real space, the order param- 
eter for this state consists of orbital currents along the bonds 
of the two dimensional square lattice of Copper atoms. Since 
the currents circulate in opposite directions in any two con- 
secutive unit cells of the lattice, the total orbital current av- 
erages to zero, and the macroscopic TR symmetry remains 
unbroken. Recently, it was shown that the admixture of a 
small d X y component to the order parameter of the DDW 
state breaks the global TR symmetry, producing a non-zero 
Kerr signal 6 in conformity with the experiments 2 . The chi- 



ral d xy + id x 2_ y 2 (d + id) density-wave state, as also the 
regular DDW and the spin density wave state, has hole and 
electron pockets as Fermi surfaces in its excitation spectra. 
Such reconstructed small Fermi pockets are consistent with 
the recently observed quantum oscillation in high magnetic 
fields in underdoped YBCO 7,8,9 ' 10,11 . In this paper, we discuss 
an intrinsic anomalous Nernst effect induced by the d + id 
density-wave state as a direct consequence of the macroscopic 
TR symmetry breaking and the presence of the Fermi pock- 
ets. Because of the broken TR symmetry, the ordered state 
acquires a Berry curvature 12 which is sizable on the Fermi 
surfaces. It is known that the non-zero Berry curvature can 
produce the anomalous Hal l 13,14 ' 15,16 and Nerns t 17 ' 18 effects in 
ferromagnets. We focus here on the anomalous Nernst effect 
for the high-T c cuprates, because the corresponding coeffi- 
cient has been extensively measure d 19,20 ' 21 . 



Nernst signal for unconventional density waves, such as 
the DDW state, was studied earlier in Ref. [22] and refer- 
ences therein. In these papers, however, the order parameter 
of the bare DDW state was used, for which the Nernst ef- 
fect was induced by the external magnetic field. On the other 
hand, in the present paper we consider the superposition of 
two different d-wave order parameters (motivated by the PKE 
measurements^ 6 -), and the spontaneous breakdown of time re- 
versal symmetry leads to the Berry curvature, which acts as a 
magnetic field. Estimating the degree of TR symmetry break- 
ing from the PKE measurements of Refi 2 -, we calculate the 
expected anomalous Nernst signal in the underdoped phase 
of YBCO near T*. We stress that even though we model 
the pseudogap by a chiral DDW state, the basic conclusions 
are more robust; the broken TR symmetry and well-defined 
Fermi surfaces, both of which have now been experimen- 
tally verified, necessarily imply the anomalous Nernst effect 
which should be observable. Note that, recent neutron scatter- 
ing experiments 23,24 have appeared to indicate a TR breaking 
state without translational symmetry breaking in the pseudo- 
gap regime 25 . We expect an anomalous Nernst effect for such 
a state as well, if it breaks the TR symmetry globally. 



II. BERRY CURVATURE OF THE CHIRAL DDW STATE 



The order parameter of the d xy + id x 2_ y 2 density wave 
stated is a combination of two density waves with different 



angular patterns 



k+Qa C k/5) = ( A k + iW k ) S a/3 , 



(1) 



where c\ c are the electron creation and annihilation opera- 
tors on the 2D square lattice of Copper atoms, k is a 2D mo- 
mentum, Q is the momentum space modulation vector (it, tt), 
and a, (3 are the spin indices. Wk = ^r(cosk x — cosk y ) 
and Ak = — Ao sin k x sin k y are the order parameter ampli- 
tudes of the id x 2_ y 2 and d xy density wave components, re- 
spectively. The imaginary part, iWk, of the order parameter 
breaks the microscopic TR symmetry giving rise to sponta- 
neous currents along the nearest neighbor bonds of the square 
lattice. The spontaneous currents produce a staggered mag- 
netic flux, which averages to zero on the macroscopic scale. 
The d xy component of the density wave, Ak, leads to the 
staggered modulation of the diagonal electron tunneling be- 
tween the next-nearest neighbor lattice sites. Such staggered 
modulation breaks the symmetry between the plaquettes with 
positive and negative circulation and, thus, breaks the macro- 
scopic TR symmetry. Such macroscopic TR symmetry break- 
ing may account for the nonzero PKE observed in the recent 
experiments 2 ^. 

The Hartree-Fock Hamiltonian appropriate for the mean- 
field d + id density wave is given by 



H 



E 



keRBZ 



£k ~ M D k exp {id k) 
£> k exp (-z6» k ) £k+Q - /i 



(2) 



where 



, £k is the free electron band struc- 

4i' cos k x cos k„, and 



k ^k+Q 

ture, £k = — 2t(cos k x + cos k y ) + 4f cos k x cos k y , ana /i is 
the chemical potential. The order parameter has been rewrit 
ten as Dk exp [idk] with the amplitude Dk = yH^ 
and the phase 6>k = ?r0 (— Ak) + arctan (Wk/Ak), where 
(x) is the step function. In writing the Hamiltonian, the 
first Brillouin zone has been folded to the magnetic or reduced 
Brillouin zone (RBZ) to treat the Q = (tt, tt) modulation ef- 
fectively. The energy spectrum of the Hamiltonian (ffjl con- 
tains two bands with eigenenergies E± (k) = wq ± w (k), 
where w Q (k) = - A t + (e k + £k+ Q ) 1% w (k) = y/F£ + D 2 
with F k = (e k - £k+Q) /2. 

Berry phase is a geometric phase acquired by the wave- 
function when the Hamiltonian of a physical system under- 
goes transformation along a closed contour in the parameter 
space 12 . For the d + id Hamiltonian (O, the relevant pa- 
rameter space is the space of the crystal momentum k. The 
eigenfunctions of the Hamiltonian are therefore k-dependent, 
and the overlap of two wavefunctions infinitesimally sepa- 
rated in the k-space defines the Berry phase connection Ak = 
($,Jj (k)| iVk |$n (k)}, where $„ (k) is the periodic ampli- 
tude of the Block wave function, and n is the band index. The 
Berry phase connection corresponds to an effective vector po- 
tential in the momentum space, and its line integration around 




FIG. 1 : (Color online) Logarithm of the Berry curvature $7+ (k) plot- 
ted on the Brillouin zone. The Berry curvature is sharply peaked at 
the points (±-§, ±f )• The ellipses and the half-circles are the hole 
and the electron pockets of the d + id state, respectively, t = 0.3 eV, 
t' = 0.09 eV, /i = -0.26 eV, W = 0.08 eV, A = 0.004 eV. 



a close path gives the Berry phase. The Berry curvature, the 
Berry phase per unit area in the k space, is given by 

n n (k) = V k xA k , A k = ($t(k)|iV k |$„(k)>. (3) 

The Berry curvature Jl„ (k), thus, acts as an effective mag- 
netic field in the momentum space and enters in the equa- 
tions of motion of the wavepacket. For a system invariant 
under both time reversal and spatial inversion symmetries, 
the Berry curvature l~i„ (k) = for every k. However, the 
d + id density-wave state breaks the macroscopic TR symme- 
try, therefore the Berry curvature can acquire non-zero values. 

To calculate the non-zero Berry curvature, we find the 
eigenstates of the Hamiltonian (O, which are given by 
$±(k) = (u ± {k)e l6k / 2 ,v ± (k)e- l9k / 2 ), where + and 
— correspond to the upper and lower bands with the en- 
ergy dispersion E + (k) and E_ (k), respectively. The co- 
efficients u± (k) and v± (k) in the eigenstates $± (k) are 
straightforwardly obtained from the matrix (f2j). Substitut- 
ing the eigenstates <£>± (k) into Eq. (O, we find Q± (k) = 
-iV k (u± (k) - v 2 ± (k)) x Vk#k- In the pure DDW state, 
#k = 7r/2 is a constant, therefore Q± (k) = and there are 
no Berry phase effects. However, in the d + id density-wave 
state, the phase 9k = irQ (— Ak) + arctan ( Wk/Ak) depends 
on the values of order parameters Wk and A&, and can vary 
in the k space, therefore £l± (k) can acquire nonzero values. 
Because the momentum k is restricted to the xy plane, only 
the z component of fl± (k) can be nonzero, which we will de- 
note as fl± (k). After some straightforward algebra, we find 



Ot (k) 



1 



2w 3 (k) 
w 3 (k) 



<9w fc dw k 

x 



dk 7 



dk„ 



(4) 



(sin 2 k y + cos 2 k y sin 2 k x ) (5) 



where Wk is a three-component vector, Wk = 
(— Ak, — Wk, Fk), and it enters into the Hamiltonian 
density in (f2]i as H = wq! + Wk • Ti. Here fi (i — I, 2, 3) 



3 



are the Pauli matrices and / is the 2x2 unit matrix operating 
on the spinors We see from Eq. © that the Berry 

curvature is nonzero only when the amplitudes Ao and Wq 
of the d xy and id x 2_ y 2 order parameters are both nonzero. 
The Berry curvatures have opposite signs in the upper and 
the lower bands: £l+(k) = — fi-(k). In Fig. Q] we plot 
the Berry curvature 17+ with respect to the momentum k 
for a set of parameters in the d + id state. We see that S1+ 
peaks at (±f , ±f ), where w (k) reaches the minimum and 
the corresponding points in the k space are the points of 
near degeneracy between the two bands. The value of il+ 
decreases dramatically along slim ellipses whose long axes 
lay on the RBZ boundary lines k y ± k x = ±tt where w (k) 
and the band splitting are the smallest. The peaks of the 
Berry curvature il± correspond to magnetic monopoles in the 
momentum spaced. 



III. ANOMALOUS NERNST EFFECT IN THE CHIRAL 
DDW STATE 
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FIG. 2: (Color online) The contour plot of the integrand of Eq. ((9} 
for the anomalous Nernst conductivity. Fhe plotting parameters are 
the same as those in Fig. 1 except A = 0.0008 eV. T = 130 K. 
Fhe main contribution to the anomalous Nemst signal comes from 
the lower band (hole pockets). Fhe upper band (electron pockets) 
gives a negligible contribution. 



In the experiments to observe Nernst effec t 17 ' 20 ' 21 , a tem- 
perature gradient — VT, applied along, say, the x direction 
produces a measurable transverse electric field. Fhe charge 
current along x driven by — VT is balanced by a backflow 
current produced by an electric field E. Fhe total charge 
current in the presence of E and — VT is thus given by, 
•h = a ijEj + OLij (—djT), where <xy and ay are the elec- 
tric and the thermoelectric conductivity tensors, respectively. 
In the experiments, J is set to zero and the Nernst signal, de- 
fined as 



e N = Ey/ | VT| = pa xy - S tan9 H , 



(6) 



is measured, where a xy is Nernst conductivity defined via the 
relation J x = a xy (—d y T) in the absence of the electric field, 
p = l/u xx is the longitudinal resistance, S = E x j |VT| = 
pa xx is the thermopower, and tan## = a xy /a xx is the Hall 
angle. For a relatively modest hole concentration away from 
the severely underdoped regime in the cuprates, the second 
term in Eq. (O is experimentally observed to be small^P.. As 
long as the second term is small, pa xy completely defines the 
Nernst signal, but in the most general case one should extract 
a xy from the experimental data, as it was done in Refs i 17 ' 20 , 
to compare with our theory. 

Fhe Berry-phase effects have found much success in 
explaining the anomalous Hall and Nernst effects in 
ferromagnet s 1 3 ' 1 6 ' 17 ' 1 8 . In the presence of an external elec- 
tric field E along the x direction, the anomalous Hall current 
is along the transverse y direction. Fhe anomalous DC Hall 
conductivity is found to be 



>:ry 



RBZ 



dk^ dky 

(27T) 2 



[/(£_ (ft))-/ (25+ (*))], 



(7) 

where / (E n ) — 1/ (1 + exp {(3E n )) is the Fermi distribution 
function at a temperature T, (3 — 1/fteT and we have used 
S1+ = — f2_. Eq. (Q agrees with the DC Hall conductivity 
of the d + id density wave obtained earlier using a different 



approac h 6,27 ' 28 . For half filling (t' = 0, [i = 0), when the 
system is a band insulator, its value is quantized, e 2 /27r? i 27,28 
per spin component. It changes continuously as the system 
deviates from half filling and the Fermi pockets appear 6 . We 
will see below that the anomalous Nernst effect is zero in the 
case of half filling and becomes non-zero only when there are 
hole and electron pockets in the spectrum. 

In order to obtain the coefficient a xy , it is more conve- 
nient to calculate the coeffecient a xy , which determines the 
transverse heat current J' 1 in response to the electric field E: 
J x = a xy E y . It is related to a. xy by the Onsager relation 
a xy = Ta xy - l& ' 29 . In the presence of the Berry curvature and 
the electric field, the electron velocity acquires the additional 
anomalous term fiv^ = eE x f2(k) 17 ' 18 . Multiplying this ve- 
locity by the entropy density of the electron gas, we obtain the 
coefficient for the transverse heat current: 



°«- Ta *y- pn^L 



n= ±JKBZ (2?r) 



ft n (k)s„(k). (8) 



Here s(k) = -/ k ln/ k - (1 - /k) ln(l - / k ) is the entropy 
density of the electron gas, /k = f[E n (k)] is the Fermi dis- 
tribution function, and the sum is taken over both bands. Us- 
ing the explicit expression for the Fermi distribution function, 
Eq. ((HJ can be transformed to the following form 



€ 1 v ^ 



dkj> dky 



f2 n x 



(9) 



{E n (k) / (E n (k)) - fc B Tl0g [1 - / (E n (k))]} 



Eq. (O coincides with the corresponding expression derived in 
Ref3 using the semiclassical wavepacket methods and taking 
into accound the orbital magnetization of the carriers 3 -. Re- 
lation of the transverse heat current to the entropy flow was 
also discussed in Refs^ anctf . At T = 0, the carrier entropy 
is zero, s„(k) = 0, so there is no heat current and a xy = 0. 
At T ^ 0, we first consider the simple case with t' = and 



4 



10.6 




FIG. 3: (Color online) Plot of the Nernst signal, see Eq. l|6}, ver- 
sus temperature T. The unit on the y-axis is nV/K. The inset, 
with temperature independent order parameters Wo = 0.08 eV, 
Ao = 0.1% Wo eV, shows that the Nernst signal satisfies the Mott 
relation dlOt . In the inset, dotted line corresponds to the Mott relation 
l |10| > and the open circles are obtained through numerical integration 
of Eq. {9j. 

fi = 0, that is, the lower and the upper bands are symmetric 
with E + (k) = — E- (k). It is easy to check that in this case 
a xy = because £1 + = — f2_. 

In the general case, the entropy s n (k) has peaks on the 
Fermi surface and decreases dramatically away from the 
Fermi surface. Because the Berry curvature f2± peaks along 
the RBZ boundary k x ± k y = ±7r, the integrand in Eq. (0 
peaks at the intersections of the Fermi surface and the RBZ 
boundary, the so called "hot spots''- 3 -^, which were shown 
earlier to be important in the calculations of the Hall coeffi- 
cient in the DDW state 33 . These peaks are clearly seen in Fig. 

El 

From Eqs. 0, (8), and (0, we can show that, at low 
temperatures, the Nernst conductivity a xy is related to the 
zero temperature Hall conductivity a xy through the Mott 
relation 34 , which yields 

= f^da^ T 
Se dfj, 

Here the derivative of a xy leads to = 

-T Jrbz [S {E - ) ~ S {E+)] - Here " HE±) 

are the delta functions. Therefore the integrand is nonzero 
only at the boundary lines of the hole and electron pockets. 
In the case of a band insulator (t' = 0, fj, = 0) that does not 
contain the Fermi pockets, a' (//) = and the anomalous 
Nernst conductivity a xy = 0, even though the DC Hall 



conductivity, Eq. 0, is non-zero and, in fact, is quantized 6 . 
This is because the quantum Hall current carries no entropy. 

For crude estimate of the Nernst signal, we choose a set of 
parameters appropriate for the underdoped YBCO 35 , t = 0.3 
eV, t' = 0.09 eV, p = -0.26 eV (corresponding to the 
hole doping of about 10%), d = 1.17 nm (the distance be- 
tween consecutive 2D layers), p ~ 3 mffcm, Wo (T) = 

0.1 (1 - T/T^) 1/2 eV, A (T) = 0.0001 (1 - T/T£) 1/2 eV, 
and numerically integrate Eq. (0, where we made reasonable 
assumptions about the transition temperatures, Tw w 150 
K 2 and T£ ~ 250 K. In Fig. we plot pa xy in a tem- 
perature regime that is below T^, but much higher than the 
superconducting transition temperature T c w 80 K 2 -. Here 
we have multiplied the results by 2 to account for the con- 
tributions from two spin components. As temperature drops 
from T^r, the order parameter Wo (T) grows, leading to the 
increase of the Nernst signal. Close to T c , the Nernst effect 
would be dominated by the mobile vortices and our calcu- 
lations do not apply there. The estimated value of pa xy at 
T ~ 130 K is about 10 nV/K. This value is about 10% of the 
experimentally observed Nernst signals in underdoped LSCO 
and BSCCO 21 at temperatures much higher than the super- 
conducting T c . Note that the spontaneous Nernst signal dis- 
cussed above may not be observable through the DC current 
measurements 21 without a non-zero magnetic field because of 
the macroscopic domains with opposite chiralities present in 
a sample at the zero magnetic field. 



IV. CONCLUSION 

In summary, we discuss the non-zero Berry curvature in 
the d + id density-wave state, which was proposed earlier 5 
to explain the time-reversal symmetry breaking 2 in the pseu- 
dogap phase of the high T c superconductor YBCO. We show 
that the nonzero Berry curvature, arising out of the broken 
time-reversal invariance, and the existence of Fermi pockets in 
the cuprates directly imply an anomalous Nernst effect which 
should be measurable. We note that measurable Nernst signals 
have been found in underdoped LSCO and BSCCO 21 even at 
temperatures much higher than T c , and we propose that a TRS 
breaking state, such as the chiral DDW state, may be the origin 
of these signals. The anomalous Nernst effect at the pseudo- 
gap temperatures will constitute a further proof of an ordered 
state, with broken time-reversal invariance, to be responsible 
for the pseudogap phenomena in the cuprates. 

This work is supported by ARO-DARPA and LPS-CMTC. 
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